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ABSTRACT 

We analyze high quahty, complete stellar catalogs for four young (roughly 1 Myr) and nearby (within 
~300 pc) star-forming regions: Taurus, LupusS, Chal, and IC348, which have been previously shown 
to have stellar groups whose properties are similar to those of larger clusters such as the ONC. We 
find that stars at higher stellar surface densities within a region or belonging to groups tend to have 
a relative excess of more massive stars, over a wide range of masses. We find statistically significant 
evidence for this result in Taurus and IC348 as well as the ONC. These differences correspond to 
having typically a ~ 10 - 20% higher mean mass in the more clustered environment. Stars in Chal 
show no evidence for a trend with either surface density or grouped status, and there are too few stars 
in LupusB to make any definitive interpretation. Models of clustered star formation do not typically 
extend to sufficiently low masses or small group sizes in order for their predictions to be tested but 
our results suggest that this regime is important to consider. 
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1. INTRODUCTION 

Does the distribution of masses of stars forming in 
isolation differ from those forming in clustered environ- 
ments? The initial mass function inferred from local 
field stars appear s to be consistent with th at seen in 
clusters (e.g.. .Bastian. Covev. fc Meveiil2010[ ). however, 
the local field star population is composed of a combi- 
nation of stars which formed in isolation, stars which 
dispersed from unbound clusters, and stars which were 
ejected from bound clusters. Most st ars are b elieved to 
form in clustered environments (e.g., iLada et al 2006), 
therefore stars which formed in isolation or small groups 
may not be the dominant contributor to the local field 
star population. 

Differences in the distribution of masses of stars form- 
ing in clusters versus isolation may be expected, par- 
ticularly at higher masses. Massive stars are known to 
form at least primarily within clusters, and it is uncer- 
tain both observationally and theoretically whether mas- 
sive stars can ever form aloncQ. In the competitive ac- 
cretion scenario, the most massive stars start to form 
early and spend most of their evolution in high density 
environments in order to accrete sufficient mass to be- 
come a massive star; lower mass st ars later form in th e 
gas around the massive stars (e.g. iBonnell et al.l 120081 ). 
Any massive stars found in isolation would therefore have 
moved there after accretion or have had their low-mass 
companions dispersed, rather than having formed in iso- 
lation. Under the monolithic collapse scenario, stars 
form from the frag mentation of their natal core (e.g. 
iMcKee fc TanI [20031 ). Depending on the physical condi- 
tions, a massive core could fragment to produce a massive 
star and / or a group of lower mass stars. Massive iso- 
lated stars are therefore not explicitly prohibited from 
forming by the model, although they may be unlikely. 
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^ The term "massive" typically refers to an OB star, with a mass 
of several Mq, however, in this paper, we also consider slightly 
lower mass stars, with spectral types as late as G. 



Cores with column densities i n excess of 1 g cm~^, where 
fragmentation is suppressed (jMcKee fc Tan|[2003[ ). have 
only been observed in highly clu stered environ ments. In 
the stationary accretion model of lMversI ()2009b[ ). massive 
stars accrete a significant fraction of their mass from the 
'clump' material beyond their own natal core. To accrete 
sufficient mass within a sufficiently short time, the clump 
gas must be denser than in isolated regions. Low-mass 
stars, on the other hand, could easily form in the sur- 
rounding lower density, filamentary distributions of gas; 
it would be difficult to form a massive star in isolation. 

Observationally, few, if any, massive stars appear to 
have formed in isolation. An upper limit of isolated mas- 
sive O stars in our galaxy was measured bv lde Wit et all 
(|200l 12001 to be 4± 2%, although this search was lim- 
ited to massive (O and B star) companions. The appar- 
ently isolated O stars may in fact belong to small clusters 
wi th lower mass comTjanions , as suggested by the models 
of IParker fc GoodwinI (l2007f ). 

Clustered star formation a ppears to extend do wn to 
surprisingly small size scales. iKirk fc MversI (|201lL here- 
after Paper I) analyzed fourteen small stellar groupings, 
typically with 20-40 members, in four young, nearby star- 
forming regions where deep spectroscopic catalogs are 
available. Despite possessing surface densities an order 
of magnitude lower than the standard cluster, iPaper II 
showed that these small stellar groups shared many of 
the properties associated with clusters. These properties 
include a correlation of the mass of the maximum mass 
member with the total group mass, and the central loca- 
tion of the most massive member. The unique advantage 
of studying clustered star formation in these small stellar 
groups is the lack of source confusion (due to a combi- 
nation of the lower source densities and closer distances 
compared to typical clusters), the inclusion of very low 
mass members (the catalogs are complete to late M spec- 
tral type), and the young age (roughly 1 Myr). Since 
only approximately half of the stars in each region are 
identified as belonging to stellar groups, while the other 
half of the stars are found in relatively isolated environ- 
ments, these data also provide an excellent opportunity 
to study differences between the clustered and isolated 
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young stellar populations, in particular the distribution 
of masses. 

Therefore, in this paper we investigate the mass dis- 
tributions of clustered and isolated populations of young 
stars. We address the question of whether more mas- 
sive stars preferentially form in clustered environments. 
Our main conclusion is that higher mass stars are pref- 
erentially found in clustered environments, over a wide 
range of stellar masses. In order to minimize uncertain- 
ties introduced by estimating stellar masses, much of our 
analysis is performed comparing the spectral types of the 
stars, which have been accurately measured. 

In Section 2, we describe the stellar catalogs and the 
completeness levels in each. In Sections 3 and 4, we 
compare the distributions of spectral types for stars in 
high and low surface density environments (Section 3) 
and those belonging to groups versus isolated stars (Sec- 
tion 4). We make comparisons to model predictions and 
various observations in Section 5, discuss the results and 
interpretations in Section 6, and conclude in Section 7. 
Stellar motion is analyzed in Appendix [Xj and several 
additional statistical tests are discussed in Appendix |B] 

2. DATA 

Currently, four nearby, young star-forming regions 
(Taurus, Lupus3, Chal, and IC348) have excellent stel- 
lar catalogues with completeness to around 0.02 Mq or 
type M 8. 5. The Tau rus data was originally compiled in 
iLuhman et all (|201Q ). while the Lupus 3 data is a combi- 
nation of several source lists given in iComeronI (|2008f ). 
In our analysis here, we exclude the six stars which lie 
at a right ascension of less than 16h 05m and a declina- 
tion below 39° 50', as we expect the stellar population is 
incomplete he re - t his area was covered by only one the 
source lists in lCome ron (20081) f h is Tab le 8). The Chal 
data is primarily from iLuhmanI (|2007| ) , and the IC348 
data is primarily a co mbination of the cat alogs given in 
iLada et all (|2006D and lMuench et all (|2007[ ). In all cases, 
the spectral types of the stars were determined spectro- 
scopically, with a typical uncertainty of roughly half a 
spectal subtype. The total number of stars in each region 
above the spectral completeness limit is 344 in Taurus, 
58 in Lupus3, 215 in Chal, and 349 in IC348. These dat a 
are discussed in more detail in Appendix A of iPaper^ . 

In lPaper^ , we identified stellar groups based on a min- 
imal spanning tree (MST) analysis. In a MST, all stars 
are linked together by their nearest neighbours in a tree 
diagram; groups are defined where all members are con- 
nected by separations less than th e critical length, Lcrit- 
Following iGutermuth et all (|2009D . Lcrit was determined 
based on th e distrib ution of all branch lengths. Ap- 
pendix D of iPapem examines the effect of Lcrit on the 
groups identified, and shows that group membership is 
little affected by var iations o f Lcrit within the range of 
errors expected. In iPaper J . we set a minimum group 
size of more than ten members, based on a visual exam- 
ination of the MST g roup ings. Unlike the 0-star studies 
of lde Wit etlH (|2004f ) and lLamb et~ail (|2010f ). foreground 
and background contamination of these catalogues is not 
a concern, allowing groups to be identified without re- 
lying on statistical comparisons with background source 
density counts. 

Masses were estimated based on a combination of 
stellar evolutionary models by iPalla fc: Stahleri (jl99^), 



iBaraffe et al (|1998D . and iChabrier et all ([20000, andas- 
suming a constant age of 1 Myr. As discussed in IPaper 1 . 
these assumptions lead to an overall uncertainty in the 
mass of each star of about 50%; the ranking of the masses 
is, however, much more certain (and is exact for any sin- 
gle stellar age assumed). Where possible in the analysis 
here, comparisons are made using spectral types rather 
than masses to minimize the uncertainties. 

2.1. Orion Nebula Cluster 

In IPaper J . we also compared the properties of the 
young stellar groups to those of a typical young clus- 
ter. For thi s purp ose, we analyzed the ONC dataset of 
IHiUenbrandl (Il997l) . There, our MST algorithm identi- 
fied one very large cluster (the central cluster where the 
Trapezium stars are found), as well as five small stellar 
groups with properties similar to those fo und in the four 
nearby star forming regions. In IPaper J . we considered 
only the large ONC cluster for comparison, as it alone 
represented typical cluster properties. 

For t he ONC cluster identification and analysis in 
IPaper J . we included only sources listed with a 70% 
or higher 'probability of membership ' (base d on proper 
motion observations) in IHiUenbrandl (|1997l ). We used 
this conservative cut to prevent potential contamination 
by non-ONC mem bers. Many of the fainter sources in 
IHiUenbrandl (jl997D . however, did not have proper motion 
data available at the time, and hence had no probability 
of membership measure. IHiUenbrandl (|1997| ) estimated 
that the majority of sources without proper motion data 
were likely bona fide cluster members, based on the frac- 
tion of known members as a function of both the sepa- 
ration from the cluster centre and spectral type. 

A comparison of the properties of isolated and clus- 
tered stars could be strongly sensitive to variable com- 
pleteness, both spatially and spectrally. In order to avoid 
any potential bias in our results for the ONC, we there- 
fore run our analysis on two versions of the ONC catalog. 
The con servative proper motion-cut catalog analyzed in 
IPaper II will be referred to as ONC-c; additionally, we an- 
alyze the full ONC catalog with no cuts applied, which 
will be denoted as ONC-f. The former should have mini- 
mal contamination from non-ONC stars, but an irregular 
completeness at later spectral types, while the latter cat- 
alog has better completeness but a higher likelihood of 
contamination. We ran the same group-identification al- 
gorithm on both ONC catalogs separately; due to the 
higher surface density, ONC-f has a smaller Lcrit than 
ONC-c: Lc«t(ONC-f) is 0.08 pc and Lc„t(ONC-c) is 
0. 13 PC. 

IHiUenbrandl (|1997D estimated that their optical catalog 
contained roughly half the total number of cluster mem- 
bers (the others being detected only in the infrared due 
to high extinction) , and that the optical sample was rep- 
resentative of the full distribution both in terms of spatial 
distribution and spectral types. Sixty percent of the op- 
tical sample had spectral types measured, with a roughly 
uniform completeness level both spatially and as a func- 
tion of I and V band photometry (assumed to roughly 
correlate with spectral types). While the s pectral type 
determinations are therefore not complete, IHiUenbrandl 
(IT997I) argue that they should be representative of the 
full population in every respect. The ONC survey com- 
pleteness is roughly 14.5 mag in K-band and 17.5 mag in I 
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band in the least-sensitive part of the survey iHillenbrandl 
(|T997), which corresponds to roughly a mass of 0.04 to 
0.055 Mq at an age of 1 Myr or 0.06 to 0.075 at 
5 Myr. We adopt a conservative estimate of a spectral 
completeness limit of M6.5. 

In general, we find similar results using either the 
ONC-c or ONC-f catalog, suggesting the completeness 
is reasonably consistent in the ONC-c catalog. 

3. DISTRIBUTION OF SPECTRAL TYPES IN HIGH 
VERSUS LOW SURFACE DENSITY ENVIRONMENTS 

3.1. Calculating the Local Stellar Surface Density 

In order to compare the properties of stars inhabiting 
isolated versus clustered environments, a scheme to clas- 
sify each is required. One simple method which does not 
rely on the definition of stellar groups is to use the local 
stellar surface density, S. If the separation from the star 
to its nth nearest neighbour (where n is 1 for the star 
itself) is r„, then the local surface density is 



As discussed ii i iGutermuth et all ()2009f ) and 
ICasertano fc HutI (|1985D . the fractional uncertainty 
in S varies as (n — 2)°-^; higher values of n give a lower 
spatial resolu t ion, but smaller fractional uncertainty. 
iBressert et all (|2010( ) recently used this surface density 
measurement to argue that there is no distinct scale for 
YSO clustering within nearby star-forming regions. 

We calculate S for each star using both n = 4 and 9, to 
allow us to examine the dependence on n. Since [Paper ]| 
found many examples of stellar groups with roughly ten 
to twenty members, n values larger than about 10 would 
have poor sensitivity to this clustering, while n values of 
two or three could bias surface density measures for close 
visual pairs or binary stars. 

3.2. Comparison of Spectral Type Distributions 

We examine whether there is a global preference for 
earlier spectral types at higher S. Ordering all YSOs 
within each region by E (for n = 4 and 9), we compare 
the spectral type distributions for the upper and lower 
thirds of the population. Figures |T] to [3] show a compari- 
son of the cumulative spectral type distributions for the 
n = 4 case for each of the four regions, as well as the ONC 
for comparison. Each plot shows the cumulative fraction 
in both linear space (main plot) and log space (inset), 
to highlight differences as the later and earlier spectral 
types respectively. Comparisons are only made for stars 
above the completeness levels discussed in Section 2 and 
with known spectral types, as are all subsequent mea- 
sures in this paper. The error bars in the figures denote 
the Poisson \/N error for each point. Three of the re- 
gions - Taurus, IC348, and the ONC - show striking dif- 
ferences between the high and low surface density stars. 
These regions show an overabundance of stars at earlier 
spectral types in the high surface density environment 
versus the low surface density environment. This is not 
only the case at the earliest spectral types, but extends 
throughout nearly the entire range of spectral types in 
the sample. Statistical tests described below confirm this 
visual impression. The results of the statistical tests are 



given in Table 1. Similar trends, but with poorer sta- 
tistical significance were found comparing the upper and 
lower E halves of the population, likely because stars at 
intermediate values of E have similar properties, but are 
distributed in both bins. 

3.3. Statistics 

We ran three statistical tests on the distributions of 
spectral types for the grouped and isolated stars to quan- 
tify the differences visually suggested in Figures [T] to O 
The first two of these tests examine global differences 
between the stars in high and low surface density en- 
vironments, which will therefore be weighted where the 
bulk of the stars are, i.e., the later spectral types. The fi- 
nal test focusses on the differences solely at early spectral 
types. 

3.3.1. Kolmogorov-Smirnov test 

We first ran a two-sample Kolmogorov-Smirnov (KS2) 
test; the KS2 test provides an effective way to mea- 
sure whether two datasets are statistically similar, i.e., 
whether it is likely that both are drawn from the same 
distribution. The KS2 test is sensitiv e to d ifferences 
in the datasets' medians and variances ()Cono vcr 1999). 
Table 1 shows the KS2 probabilities that the high and 
low surface density stars share a common distribution 
of spectral types. The KS2 test suggests that the high 
and low surface density stars are distinct at the 95% or 
higher confidence level with either n = 4 or 9 in IC348 
and in the ONC, while in Lupus3, the populations are 
consistent with being drawn from the same distribution 
at the 96% confidence level, likely due to the small num- 
ber of stars in that region. Chal has intermediate KS2 
probabilties and so does not fall definitively into one or 
the other category. Taurus shows a strong difference in 
populations for n — 4 (95% confidence level), but not for 
n = 9, suggesting that very small stellar groupings are 
significant in the spatial distribution of stars in Taurus. 

We ran similar comparisons between stars in high and 
low surface density environments with more severe (ear- 
lier) completeness levels assumed, and find similar KS2 
values for completeness levels of 3-5 spectral subtypes 
higher in nearly every region, for both the n = 4 and 
9 surface density measures. The one exception is the 
ONC-c catalog, where KS2 values remained low only for 
a completeness level of up to two spectral subtypes ear- 
lier than what we assumed. 

The KS2 test does not have a formal mechanism to 
include measurement errors in the calculation. In order 
to assess the effect of uncertainty in the spectral types 
on the KS2 test statistic, we therefore ran a series of tri- 
als with random errors added to the measured spectral 
types. For each entry in Table 1, we ran 10,000 trials 
with added spectral type uncertainties of half a subtype, 
and calculated the resulting KS2 test statistic. This is 
summarized in Table 2, and demonstrates that the signif- 
icant KS2 test statistics are robust to uncertainty at the 
expected level. Systematic biases in the spectral types 
measured are not expected as a function of spatial clus- 
tering or spectral type (K. Luhman, private communica- 
tion). 

3.3.2. Mann-Whitney 
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While the KS2 test is an excehent way to determine 
whether two distributions are different, it does not pro- 
vide a measure for how they differ, e.g., if one of the 
distributions tends to have extra high- or low- valued 
members. To address this, we ran a second statistical 
test, the Man n- Whitney (MW ) test, also known as the 
Wilcoxon test ()Conoveiill999[ ). In this test, two datasets 
are compared based on their relative ranking (i.e., the 
ordering of the values of both datasets combined); tied 
values are each assigned the average rank of the ties. The 
basic premise of the test is that two sets of data drawn 
from the same overall distribution will have roughly equal 
total ranks, whereas if one data set tends to have higher 
values, then the total of its ranks will tend to be lower 
(since a rank of 1 denotes the largest value). The MW 
test probability that the high surface density stars tend 
to have earlier spectral types than the low surface density 
stars is given in Table 1 for each region. We emphasize 
that the MW test is most sensitive to differences in the 
range of spectral types where the bulk of the stellar pop- 
ulations are, i.e., earlier spectral types refers primarily 
to those at the late-type end, around M-type. Taurus, 
IC348, and the ONC show high probabilities that the 
stars in high surface density environments tend to have 
earlier spectral types. Both LupusS and Chal have no 
strong evidence that either the higher or lower surface 
density stars tend to be of earlier spectral type. 

Where differences are seen in the spectral type distri- 
butions, we can make a rough estimate of this typical 
difference. We apply a global shift to the spectral types 
of the low surface density stars, and measure how large an 
offset may be added so that the MW test shows negligible 
statstical significance that the high surface density stars 
tend to have earlier spectral types. We find that required 
shifts in the spectral types are quite small - shifting all 
of the low surface density stars earlier by between half 
and one spectral subtype is sufficient for the MW test 
to return probabilities of 859^ or smaller that Taurus, 
IC348, and ONC stars have earlier spectal types in high 
surface density environments. 

As with the KS2 test, we ran the same statistical com- 
parisons varying the assumed completeness level of the 
sample, and found the same result, i.e., the conclusions 
drawn from the MW test remain unchanged for complete- 
ness levels up to at least two spectral sub-types earlier 
than assumed. 

Uncertainties in the spectral types measured also do 
not have an impact on the MW test statistic - we ran 
a similar test to that described in the previous section 
for the KS2 test, and find little variation in the MW test 
statistic (see Table 2). 

3.3.3. Early Type Stars Counting Test 

Neither of the above tests explicitly focusses on differ- 
ences specifically in the early type star populations at 
high and low surface densities. There are an insufficient 
number of stars in each region (save the ONC) to run the 
above two tests only on the early type stars in each en- 
vironment. Instead, we must rely on a simpler statistical 
test, comparing the total number of early type stars in 

^ The value of 85% is chosen to represent a probabihty that, 
while reasonably high, is low enough to prevent any firm conclu- 
sions. 



high and low surface density environments. If the popu- 
lations were identical, roughly the same number of early 
type stars should be found in the high and low surface 
density samples, since the total number of stars in each 
is equal. To determine whether the differences are sig- 
nificant, we compare the value to the Poisson error in 
the number of early type stars in the low surface density 
sample. We then define the excess of early type stars in 
the high surface density environment as 

771 Nearly, high sd^ ^ early .Low sd 

-Nearly / . . — 

Sji^/Nearly.high s d)'' + {y/ N ^arly ,lan. ,dV 

Nearly, high sd — ^ early , low sd ^ ' 

yj2Nearly,Low sd 

which gives roughly the number of sigma of significance. 

The value of E^ariy depends on what spectral type is 
used as the early type cutoff. Since it is not obvious, a 
priori, what spectral type should be used, we performed 
the calculation for a range of cutoffs, every half a spectral 
type from late G to early B, in order to see what type(s) 
would give the highest statistical significance to the dif- 
ference. Table 1 gives the maximum value of E^ariy we 
found, and the early type cutoff spectral type or types 
that it was found for. The four nearby regions tend to 
have Eeariy valucs of < 2a with an early type cutoff of 
around GO. The ONC had much higher Eeariy values, 
from 3 to 8 (7, when a very early type cutoff, around B5, 
was used. 

3.3.4. Statistics Summary 

In summary, we find that there are statistically sig- 
nificant differences between the spectral types of stars 
found in low and high surface density environments in 
IC348 and the ONC, as well as Taurus, when small stellar 
groupings are considered (using n = 4). In these regions, 
the higher surface density stars have an overabundance 
of earlier spectral types compared to the lower surface 
density stars. This overabundance corresponds to an av- 
erage spectral type which is half to one spectral subtype 
earlier for the high surface density stars. There are too 
few stars in Lupus3 to show any statistically significant 
differences, while the Chal stars do not follow the gen- 
eral trend. Focussing only on the early type stars, none 
of the four nearby regions have a sufficient number of 
stars to show statistically significant differences in the 
populations. Significant differences are seen, howevever, 
in th e ONC, as has also bee n noted in previous studies 
(e.g.. IHillenbrand fc Hartma nn 1998). 

3.4. Combined Distribution 

Finally, we examine the spectral type distributions for 
the combination of all four regions. Combining the spec- 
tral type distributions may increase the statistical signif- 
icance of the difference between the high and low surface 
density environments, since all regions but Chal visually 
follow the same trends. Merging the individual spectral 
type distributions can be done in two ways. The first 
method is to add together the spectral types for all of 
the YSOs classified as high surface density in their indi- 
vidual regions, and similarly for the low surface density 
stars. This is shown in the top panel of Figure |31 and the 
resulting statistical measures are given in Table 1. The 
appearance is qualitatively similar when stars in Chal 
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are excluded, however, the statistical significance of the 
trends is slightly higher in this instance (these numbers 
are also given in Table 1). 

On the other hand, the data from the four regions can 
also be combined by splitting up the stars into those 
falling in the upper and lower global thirds of the sur- 
face densities measured in all four regions. This would 
be the more appropriate combination method to adopt 
if the mechanism driving the difference between environ- 
ments was dependent on the absolute surface density of 
stars. The bottom panel of Figure |4] shows the combined 
spectral type distributions split by absolute surface den- 
sity. Note that each region populates a different fraction 
of the high and low surface density cuts; Taurus YSOs 
dominate the low surface density cut, while IC348 YSOs 
dominate the high surface density cut. The difference be- 
tween the spectral type distributions using the two com- 
bination methods is striking. While both show an excess 
of early type stars in the higher surface density environ- 
ment, the size of the difference and the range of spec- 
tral types over which the difference is present are both 
smaller when the global surface density cut was applied. 
The smaller range of the excess has a significant effect on 
the KS2 and MW statistics measured, since much of the 
population is at later spectral types where the difference 
between high and low surface densities is diminished. 

The fact that combining the data from the four re- 
gions using the relative surface density cuts rather than 
a global cut shows a much larger difference between the 
two populations suggests that the mechanism responsible 
for the difference at high and low surface densities oper- 
ates in a relative sense, rather than coming into effect at 
a set absolute surface density value. It is interesting to 
note that the global surface density comparison shows a 
much stronger difference when Taurus rather than Chal 
stars are excluded from the sample, despite Taurus stars 
on their own showing a difference while Chal stars do 
not. This appears to be the case because Taurus's early 
type stars, which are in relatively higher surface density 
environments within Taurus, are classified as being in 
low surface density environments in comparison to sur- 
face densities seen in the other regions. 

4. DISTRIBUTION OF SPECTRAL TYPES IN GROUPS 
VERSUS ISOLATION 

A second way to compare the stars is to subdivide the 
stars in each region into those that w ere found to belong 
to an MST-defined stellar group in iPaper II and those 
that were not associated with a group. As discussed in 
Section 2, groups are defined as having more than ten 
members with nearest neighbour separations less than 
Lcriu stars not associated with a group ('isolated stars') 
are therefore all stars which are connected to fewer than 
ten sources with separations of Lcrit or less. Figures[5]to 
[7] show the cumulative fraction of sources earlier than a 
given spectral type for each of the four regions as well as 
the ONC for comparison, using the same plotting con- 
ventions as Figure [TJ 

Unlike the comparisons at high and low surface den- 
sities, the grouped versus isolated cumulative distribu- 
tions tend to show significant differences over more lo- 
calized spectral types. IC348 in particular, and Taurus 
to a lesser extent, appear to show an excess of later- type 
stars (mid-M) in the groups, and this also appears in the 



ONC. The main difference between the two ONC cata- 
logs is the stars in groups with spectral types around F. 
iHillenbrandl (|1997f) found that stars around type F were 
the most likely to have small proper motion membership 
probabilities; if proper motions were preferentially mea- 
sured for stars nearer the centre of the ONC, then our 
ONC-c catalog would have a deficit of F stars in groups 
versus isolation compared to the full ONC-f catalog, as 
Figure [7] shows. 

4.1. Statistics 

We ran the same statistical comparisons as discussed 
in Section 3 on the grouped and isolated stellar distribu- 
tions. These results are summarized in Table 3. 

In general, we see similar results to the statistical tests 
for the high versus low surface density stellar distribu- 
tions, although with generally slightly poorer statistical 
significance. We again verified that the spectral com- 
pleteness limits adopted do not affect the results, for 
completeness levels at least 2-3 subtypes earlier than we 
assumed. 

Here, to compare the number of the early type stars 
in groups verus isolation, we first correct for the total 
number of stars in each category. The number of early 
type stars expected to be found in groups, based on the 
number of early type stars found in isolation, is 

M —AT Y ^tot,9rp /o\ 

early ,grp^exp — early ^isol ^ X'J ) 

i^totAsol 

and the excess of early type stars in groups over that 
expected from the isolated population is analogous to 
equation 2, but using the full error propogation required 
from equation 3. 

Similar to the high and low surface density popula- 
tions, we find relatively small excesses of early type stars. 
The excesses are typically smaller than for the high and 
low surface density populations, since the extra step of 
normalizing the total sample sizes increases the associ- 
ated error measure. Here, the largest magnitude of ex- 
cess in Chal is actually a deficit of early type stars in 
groups, although at less than Icr significance. The ONC 
stars again show the most significant excesses of early 
type stars in the groups, but only at the < 2a level. 
We tested the effect of measurement uncertainty in the 
spectral types using the same procedure as outlined in 
Section 3.3.1. The results, listed in Table 4, show that 
the KS2 and MW test statistics are robust to the spectral 
type error. 

In summary, similar to the surface density compar- 
isons, the tests show that IC348 and the ONC have statis- 
tically significant differences between the spectral types 
of stars in groups versus isolation. In Taurus, IC348, 
and the ONC, stars in groups have a global tendency to 
be at earlier spectral types than stars in isolation, again 
with a high statistical significance. This tendency is not 
obvious when examining only the earliest spectral types 
(with the Nearly test), likely due to small number statis- 
tics. The trends seen in these three regions tend to be 
slightly statistically weaker than those seen comparing 
the high and low surface density environments, and may 
be more pronounced at localized spectral types. Lupus3 
again does not have a sufficient number of stars to see 
a statistically significant difference, and Chal does not 
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appear to follow the general trend. 

4.2. Combined Distributions 

Similar to Section 3.4, we also compare the distribution 
of stellar spectral types for stars in and outside of MST 
groups in Taurus, LupusS, Chal, and IC348 combined. 
This data is shown in Figure [51 and the corresponding 
KS2 and MW statistics are given in Table 3. As ex- 
pected, combining the data from the four regions tends 
to enhance the trends seen, as is also apparent in the sta- 
tistical measures. Excluding the stars in Chal, where no 
trend is apparent, has little effect on the qualitative ap- 
pearance of Figure[Sl but does slightly further strengthen 
the statistical measures. 

4.3. Individual Taurus and Chal Groups 

While the populations of grouped stars are dominated 
by a single clustered environment in both IC348 and the 
ONC, and only one group was identified in Lupus3, in the 
Taurus grouped population, there are comparable contri- 
butions from eight different stellar groups. Here, we ex- 
amine each group independently to determine whether 
our conclusions change when groups are treated sepa- 
rately. Similarly, we compare the three groups identified 
in Chal, noting that there, one of the groups dominates 
the total number of sources by a factor of two. 

The number of stars in each group above the spectral 
completeness limit ranges from fourteen to thirty one in 
Taurus, and from twelve to eighty three in Chal, which 
are small sample sizes to search for statistically signifi- 
cant differences. Nevertheless, we compared the KS2 and 
MW test statistics for each group individually versus all 
of the isolated stars in the region. These statistics are 
given in Table 5. Results from the Nf.ariy test are not 
shown, as none of the groups have a statistically signif- 
icant excess of early type stars at the earliest spectral 
types, due to the small number of sources in each group. 

Most of the Taurus groups have KS2 probabilities in- 
dicating the single group and isolated stars are statisti- 
cally indistinguishable, and MW probabilities that nei- 
ther population has a strong tendency for more mas- 
sive members. Surprisingly, two Taurus groups did show 
strong global differences with the isolated population, de- 
spite the small number of members. Both Taurus Groups 
4 and 5 (L1551 and L1529) have KS2 probabilities of 
less than 10% and MW probabilities of > 90%, implying 
that overall the stars in these groups tend to have earlier 
spectral types than the isolated population. The spec- 
tral type distributions of the exceptional Taurus Groups 
4 and 5 are shown in Figure |9l 

Group 5 / L1529 is a small group, with its most massive 
member only of spectral type K5, and shows a substan- 
tial excess of sources relative to the isolated population 
around type MO. Group 4 / L1551, on the other hand, is 
a somewhat larger group, with a B9 / B9.5 pair of stars, 
and shows a systemati c offset of more massive stars at 
all spectral types. In IPaper J . L1551 was one of only 
two groups which did not have a centrally-located most 
massive member - the B9 and B9.5 stars fall near the 
outskirts of the group. L1551 therefore is only an exem- 
plar for the differences seen in the spectral type distri- 
butions in groups and isolation, and not also the loca- 
tion of the most massive member. In IPaper^ . we noted 



that had the group been defined with a slightly smaller 
Lcrit, the B9/B9.5 stars would not have been considered 
group members, and the next most massive member is 
centrally- located. Interestingly enough, if the B9/B9.5 
stars were excluded from the L1551 spectral type distri- 
bution, the grouped and isolated stars would still show a 
systematic difference - the absence of the two B9/B9.5 
stars is not sufficient to make the two distributions agree. 

Conversely, we test whether the excess of earlier type 
stars is confined to only these two groups in Taurus. We 
compare the spectral type distributions for the stars in 
all groups except groups 4 and 5 and the isolated stars, 
and find a diminished, but still visible, excess of early 
type stars in groups. This excess only extends down to 
about type KO, and is sufficiently small that the KS2 and 
MW tests do not yield significant results. This appears 
to indicate that the bulk of the difference in the grouped 
stars populations arises from only a few of the groups, 
although there are hints that the difference is present at 
a lower level in at least some of the other groups. 

In Chal, each individual group appears much more sim- 
ilar to the isolated distribution; none of the groups have 
the striking visual appearance of differences that are seen 
in Figure [HI Chal Group 2 has a relatively low KS2 
probability that the grouped and isolated stars are simi- 
lar (12%); this appears to be caused by a slight deficit of 
grouped stars at mid-M spectral types versus the isolated 
distribution. 

4.4. Distribution of Spectral Types with Varying Croup 

Size 

Our analysis in Section 4.1 distinguishes groups and 
isolated stars assuming that groups must have more than 
ten members. The MST group identification scheme 
used in Paper I, however, identifies all groupings of stars 
which are separated by less than Lcrit, and therefore 
identifies groupings as small as pairs of stars. T he mini- 
mum group size is somewhat arbitrary, chosen in IPaper^ 
to allow for the measurement of mass segregation. In 
light of the fact that very small groupings of stars may be 
important, particularly in Taurus (Section 3), it is useful 
to examine various minimum group sizes. A visual exam- 
ination of the YSOs suggests that a significant fraction 
of the earliest type stars not found in our iV > 10 MST 
groups are instead found within smaller MST groupings; 
there are relatively few early type stars which do not 
belong to any size of MST grouping. 

We re-ran the statistical comparisons discussed in Sec- 
tion 4.1 on grouped and isolated stars using different 
minimum group sizes. The details are discussed in Ap- 
pendix [Bl the overall conclusion is that the definition of 
a group does not have a strong infiuence on the differ- 
ences seen between the grouped and isolated stellar spec- 
tral types. Increasing the minimum group size substan- 
tially causes some of the genuine groups to be classified as 
isolated sources, diluting the difference between the two 
populations. Decreasing the minimum group size has a 
small effect on the statistical measures until the mini- 
mum group size is so small that few stars are categorized 
as isolated and the ability to statistically distinguish the 
two populations is poor. Table 6 shows the results for 
the same statistical tests as given in Table 3 but with a 
minimum group size of > 2 and A > 15. As in Sec- 
tion 3, the stars in Taurus show the largest increase in the 



7 



KS2 probability with an increase of the minimum group 
size, suggesting that stellar clustering on small scales is 
important in this region. 

5. DISTRIBUTION OF MASSES 

Comparison between our observations and other obser- 
vations, in addition to predictions from clustered star for- 
mation models are more easily made using the distribu- 
tions of mass, rather than spectral types. As illustrated 
in Sections 3.4 and 4.2, the distribution of spectral types 
combined over the four nearby star-forming regions of- 
fers a reasonable representation of the trends seen in the 
individual regions, with a generally higher statistical sig- 
nificance due to the greater number of stars compared. 
For simplicity in mass space, we make comparisons only 
of these combined datasets. Note that for the surface 
density comparison, we use the combination displayed in 
the top panel of Figure |4l i.e., each region contributes 
an equal number of YSOs to the high and low surface 
density populations^ 

As discussed in iPaperl . our simple method of assum- 
ing a constant stellar age tends to introduce biases in 
the mass distribution. To better understand this effect, 
we made comp arisons of the mass distributions shown in 
iLuhman et aj (|2009^ for Taurus, Chal, and IC348 versus 
the mass distribution which we would derive using the 
same basis set of spectral types and a constant 1 Myr age. 
We find that the assumption of a global 1 Myr age tends 
to cause an overestimation of the number of the most 
massive stars in each region (around 2 Mq and above), 
and hence an accompanying underestimate of the number 
of slighly less massive stars; there is also a smaller preva- 
lent overestimate of the number of stars with around 
0.2 Mq, and a similar underestimate at slightly lower 
masses. Assuming a constant age of 2 Myr instead, we 
find significantly improves our mass estimations above 
~1 Mq; we could not find a single age estimate to im- 
prove the smaller bias at the lower mass end. For our 
analysis of the mass distributions we therefore revise our 
masses above 1 M© to be those predicted assuming an 
age of 2 Myr, which slightly reduces the values. We stress 
that our mass distributions should still be viewed as ap- 
proximate estimations, given our simple age assumption. 

Figure [TU] shows the cumulative mass distribution for 
the regions, overlaid with predictions from several models 
(more details below). The differences between the pop- 
ulations are less visually apparent than in Figures |4] and 
[HI the difference is entirely due to the variable stretch in 
the conversion from spectral to mass space, since a single 
conversion between spectral type and mass is assumed. 
The bumps in the mass distribution are likely due to in- 
accuracies in our assumed conversion. Note that despite 
the visual appearance, our statistical tests provide the 
same results in mass-space as spectral type-space, since 
we assume a one-to-one translation between the two. To 
quantify the typical mass difference between the two dis- 
tributions, we ran a similar MW shift test to that de- 
scribed in Section 3.3.2, finding the shift which, when 
applied to the low surface density or isolated stars, re- 
turns an MW probability of ~85%. We found that for 
both Taurus and IC348, the high surface density stars 
are typically 18% more massive than the stars in low sur- 
face density environments, while stars in groups are 4% 
and 11% more massive than their isolated counterparts 



in Taurus and IC348 respectively. 

For comparison. Figure [11] shows the corresponding 
plots for the ONC-c catalog; the ONC-f catalog appears 
similar. In the ONC, th e mass estimates are expected to 
be more accurate, since iHillenbrandl ()1997l ) allowed the 
stellar ages to vary when estimating the masses. Apply- 
ing the same MW shift test as above, we find that the 
stars in high surface density environments are roughly 8% 
more massive than the stars in low surface density envi- 
ronments, while the stars in groups are 6% more massive 
than their isolated counterparts. 

5.1. Effect of Single Age on Statistics 

The statistics discussed above were calculated assum- 
ing a constant stellar age, however, variations in the stel- 
lar ages could affect the values quoted. Here, we discuss 
additional tests to investigate this effect. 

5.1.1. KS2 and MW Test Statistics 

As noted above, with the assumption of (any) constant 
stellar age, the KS2 and MW test statistics are identical 
to those found for the spectral type distributions, since 
there is a one-to-one mapping between the two. Allow- 
ing for age variations could change the shape of the mass 
distributions. In particular, having systematically older 
stars the more clustered environments would tend to de- 
crease their masses, and therefore lessen the difference 
seen between the more and less clustered populations 
(i.e., high versus low surface density or grouped versus 
isolated). We ran tests varying the stellar ages adopted 
for the clustered stars to determine the age difference 
size necessary to make the two populations appear simi- 
lar. Assuming that the less clustered population is 2 Myr 
old, then we find that the more and less clustered popu- 
lations can match at the very high mass tail (stars above 
~2 Mq) only when the more clustered stars are ~5 Myr 
old. Even this large an age difference, however, is in- 
sufficient to decrease the excess of more massive stars 
in the more clustered population in the ~0.2 - 2 Mq. 
This remains true even if the less clustered population 
is reduced to an age of 1 Myr. Since stellar evolution 
tracks run roughly vertically for masses below of order 
1 Mq , changes in age have almost no impact on the stel- 
lar masses inferred from the spectral type. 

Since age thus primarily affects only the highest mass 
stars which are few in number, the KS2 and MW statisti- 
cal tests are largely impervious to systematic changes in 
age. For any reasonable assumption about ages, the KS2 
test statistic remains small (5% or less) in all cases where 
it was small under the assumption of a global constant 
age. Age has a slightly higher, but still small, effect on 
the MW statistic. In Taurus, the MW statistic drops by 
less than 2% for more clustered stars which are 1 Myr 
older than the less clustered stars; only when the more 
clustered stars are set to 5 Myr while the less clustered 
stars are 1 or 2 Myr does the MW statistic drop below 
90% in some cases. Even in this drastic case, the MW 
statistic for the n = 4 surface density comparison only 
changes from 99% to 98%. In IC348, the effect is slightly 
larger, due to a greater number of sources at the higher 
mass ranges which are effected by the age assumed. Still, 
the MW test statistic drops by at most a few percent 
when the more clustered stars are < 1 Myr older than 
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the less clustered stars. All of the IC348 MW test statis- 
tics remain above, and normally well above, 90% while 
the more clustered stars are younger than 5 Myr. 

5.1.2. Mass Shift Values 

We find that the mass shift values quoted are only very 
weakly dependent on the stellar ages assumed. We ran 
similar comparisons between the more and less clustered 
stars, and found that the mass shifts were roughly identi- 
cal for any constant age adopted for both popuations. If 
the more and less clustered populations have ages differ- 
ent by less than 0.5 Myr, again the mass shifts measured 
remain virtually unchanged. For more clustered stars 
which are older (/younger) by 1 Myr, the mass shifts are 
decreased (/increased) by at most 2%. 

5.1.3. Summary 

The statistical results presented for the mass distribu- 
tions are largely insensitive to variations in the age as- 
sumed, despite age being a strong factor of uncertainty 
in the absolute values of the masses estimated. In or- 
der to decrease the statistical significance of the results 
presented in this section, stars would need to be sys- 
tematically and noticeably older (a > 1 Myr difference) 
than stars in the less clustered environments. There is 
no clear physical mechanism which would cause more 
clustered stars to be systematically older by such an 
amount. In fact, clustering is often seen to be tighter fo r 
younger stars than older ones (e.g.. iLuhman et a]||2010[ ). 
Systematically younger stars in clusters would tend to 
strengthen our statistical measures and would tend to 
increase the mass shift measured between the more and 
less clustered environments. 

5.2. Comparison to Model Predictions 

Here, we compare the observed cumulative distribu- 
tion of masses to two recent model predictions - the 
first, analysis of compe tetive accretion simulations in 



iMaschberger et al (|201C , hereafter MIO] , and t he second , 
analytic predictions bv iKrumholz et all ()2010l hereafter 
KIO) based on the monolithic collapse scenario. 

5.2.1. Competitive Accretion Simulations 

iMld analyzed the clustering properties of stars (rep- 
resented by sink particles) formed in th e competi- 
tive accretion simul ations of iBonnell et al.l |2008) and 
iBonnell et al.l (|2003f ). with a focus on the evolution of 
clusters and their most massive members. Clusters were 
identified using the MST formalism, with Lcrit set to 
0.025 pc through visual inspection and a minimum size of 
thirteen members. For comparison, in Papejll, we deter- 
mined Lcrit based on the distribution of branch lengths, 
with Lcrit ranging from 0.52 pc in the more dispersed 
Taurus star s to 0.083 pc in the more tightly clustered 
IC348. The IMld critical length for highly clustered sim- 
ulations therefore seems reasonably consistent with our 
definitions; the minimum number of group members is 
al so sim ilar. 

ImToI found that about 60% of all sources form within 
clusters, and of those sources, most tended to form with 
the same radial distribution as already existed in the 
cluster. The most massive sources (defined as any source 
with a final mass greater than about 1 M©), on the 



other hand, tended not to form in clusters - rather, clus- 
ters tended to later form around a pre-existing massive 
source. Often, this started out as a small 'sub-cluster', 
which then successively merged with other small 'sub- 
clusters', and eventually formed a massive cluster, with 
the massive stars from each sub-cluster ending up in 
the centre of the final system. This evolution occurred 
rapidly, as the simulat ions were only run fo r 0.5 Myr. 

At the end of the IBonnell et al.l (|2008l ) simulation, 
the distribution of masses of sink particles found within 
groups / clusters and in isolation and the fits listed in Ta- 
ble 1 of Mil The sink particle masses are upper limits to 
the masses expected in the stars - the authors noted that 
effects such as stellar feedback and binary fragmentation 
were not included in the simulation, and would serve to 
decrease the resultant stellar mass. With that caveat. 
Figure 110! shows the sink particle mass distributio ns in 
IMIO. overlaid on the observed stellar distributions. IMld 
only measure the mass distributions down to 0.8 Mq, 
which corresponds to roughly 15% of the total number 
of source s in all of our observations; the lines plotted 
show the IMld results scaled down by this factor. 

A comparison between the scaled-down MIO distribu- 
tions and the observed ones show that the two are in rea- 
sonable agreement, given both the small number statis- 
tics (and estimated mass uncer tainty) in the observa- 
tions, as well as the errors IMld estimate in their mass 
function parameters (not plotted). It is also interesting 
to note that the IMld results show a substantial distinc- 
tion between the clustered and isolated sinks, with most 
of the massive sinks being found in clusters. We observe a 
much smaller difference between the stars in groups and 
isolation. It is possible, however, that some processes 
which can reduce the final stellar mass from that in the 
simulated sink particles (such as stellar feedback) operate 
more strongly in a clustered environme nt, w hich would 
decrease the difference between the two IMld curves. 

5.2.2. Monoltthtc Collapse 

An analysis similar to lMld has not yet been performed 
for the monolithic collapse scenario, making a direct com- 
parison to our observati ons c hallenging. The best avail- 
able predictions are from lKlOl . who construct several sim- 
ple analytic models to describe the distribution of final 
stellar masses based on results from several star forma- 
tion simulations which include radiative feedback. 

Based on both t he simulat ions and a prio r theoretical 
prediction (|Krumholz fc McKeell2008l ). [KIOI find that in 
high density environments (above ~ 1 g cm~^), fragmen- 
tation of massive dense cores is largely inhibited, allowing 
for the formation of massive stars. In lower density en- 
vironments, fragmentation occurs much more regularly, 
leading to smaller stellar masses arising from a given 
core. These r esults were used as the basis for an ana- 
lytic model in iKld which predicts the final distribution 
of stellar mass es for med in both high and low density 
environments. iKld ass ume the initial dense core mass 
distribution follows the iChabriei] (|2005f ) IMF, scaled up 
by a factor of three, and a star formation efficiency for 
any core fragment of a factor of 1/3. In the high density 
case (model 'fH'), they assume no fragmentation, i.e., 
each cor e forms one star , and the mass distribution mir- 
rors the lChabrierl (|2005f ) IMF. In the low density regime 
('fL'), several different cases of fragmentation are consid- 



9 



ered, with varying amounts of mass in each fragm ent. 

The high and low density regimes of the |K13 models 
do not have a direct corresponendence to our observa- 
tions, as the spatial distribution of high and low density 
cores was not specified in the model. While a direct 
comparison therefore cannot be made, the observed dis- 
tribution of masses for grouped and isolated stars should 
fall within the range spanned by the m odels. Figure [TOl 
shows the range spanned by the iKld models at 1 and 
10 M0. (Note that these models span all masses, and do 
not require norm alization for comparison with observa- 
tions as the lMlOl measurements do.) 

The most noticeable feature of the iKld models is that 
the fraction of stars expected at 10 Mq and above is 
relatively high for both low and high density environ- 
ments. All of the IKIOI models predict more stars at 
10 Mq and above than even the IMIOI distribution for 
groups; the iKlOi range lies well above the observed value 
for the four nearby star-forming regions, and is barely 
consistent with the high surface density / grouped stel- 
lar population (only) i n the ONC. The predicted fraction 
of massive stars in the lKld models could be decreased by 
assuming a smaller maximum core mass, or, in the low 
density case, by assuming either a smaller fraction of the 
mass forms the main source, or that the core fragments 
into a larger number of pieces. 

5.3. Comparison to Field Star MF 

Finally, we compare to observed mass functions. Mass 
functions are often presented in differential form, i.e., 
number per mass bin, so we follow that convention here. 
Figure [T^] shows the differential mass function for high 
and low surface density stars (top panel) and grouped 
and isolated stars (bottom panel) . The general tendency 
is for the high surface density stars to be more numer- 
ous than the low surface density stars at the upper mass 
end, and the reverse at the low mass end. Similarly, the 
grouped stars tend to be more numerous than the iso- 
lated stars at high mass, and less numerous at low mass. 
This common behaviour is more prominent for stars at 
high and low surface density than for stars in groups and 
isolation, as was seen earlier. 

We note that inaccuracies due to our assumption of 
a constant stellar age are more apparent in differential 
than cumulative form. In Figure llOi assuming an age 
of 1 Myr or 2 Myr when estimating the masses above 
1 M0 creates barely discernable differences in cumulative 
mass distributions. In the differential mass distribution, 
however, there is a noticeable shift of objects from the 
highest mass bin to the second highest mass bin when 
the ages are assumed to be older. Most of these stars 
which shift bins have masses little above the boundary 
between the two bins with an assumed age of 1 Myr; 
ou r bin sizes and bound aries were chosen to match those 
in iLuhman et all ()2009() to facilitate easier comparison 
between the two. 

Keeping in mind that the estimated masses are still 
somewhat uncertain, we compare the distribution of 
masses in the four nearby star-forming regions to the 
field star mass function. The field star mass function pro- 
vides an approximate measure of the current, nearby iso- 
lated star-forming population. These stars are expected 
to have formed over a range of times, and therefore at the 
present day will be somewhat depleted of their higher- 



mass members, relative to the initial mass function. Fur- 
thermore, the field star population is expected to consist 
of a combination of stars which formed in isolation, with 
those that formed in small groups or clusters that have 
since dispersed or ejected some of their members. The 
field star mass function would thus be expected to be 
formed from a combination of our more and less clus- 
tered stellar populations^ 

Figure [12] shows the iReid et all (|2002D field star mass 
function, corrected for the effects of stellar evolution 
(their Figure 14, bottom panel); this field star mass func- 
tion spans a larger ra nge in masse s than some more recent 
observations such as iCovev et all (|2008f) . The evolution- 
corrected field star mass function is consiste nt, within 
errors, with other determinations of the IMF ()Reid et all 
[2OOI 1. While the mass distribution of the stars in high 
surface density environments appea rs to b e skewed to 
slightly higher masses than the iReid et all p002, ) field 
star mass function, they are consistent within errors. The 
low surface density, grouped, and isolated stellar mass 
dis tributions all app ear to be in closer agreement with 
the IReid et all (|2002l ) field star mass function; the largest 
discrepancy is around 0.2 M©, where we earlier noted 
our mass estimates were biased relative to those where a 
range of ages are used. 

5.4. Maximum Versus Total Group Mass 

Another property of the groups examined in iPaper 11 
was the relationship between the mass of the most mas- 
sive group member and the total group mass. In large 
clusters, the mass of the maximum mass member and 
the t otal cluster mass are correlated fe.g.. iWeidner et all 
I20T0I ). At high masses, the flattening of the relation- 
ship between maximum and total mass was interpreted 
as evidence that a maximum stellar mass exists. At 
lower masses, the relationship has a constant slope which 
is roughly consistent with the relationship expected for 
cluster memb e rs ran domly sampling the IMF, although 
iWeidner et all (12010') find that the scatter in the clus- 
ter data is smaller tha n expec ted for pure random sam- 
pling. We showed in IPaper II that the relationship be- 
tween maximum and total mass continues to the smallest 
of the groups we identified. Here, we examine whether 
the trend is also present in even smaller grouping of stars, 
and whether there is a minimum group size at which the 
correlation is seen. 

The top panel of Figure |T3| shows the relationship be- 
tween maximum and total stellar mass for all group- 
ings of six or more stars in Taurus, LupusS, Chal, and 
IC348, while the bottom panel shows the same relation- 
ship for groupings of two or more stars. In both pan- 
els i_th£_grej;_df^nion^ show the cluster data compiled 
bv IWeidner et all ()2010f ). while the coloured circles show 
the stellar groupings we identified with the MST. The 
dotted line shows t he approximate relationship fit by 
IWeidner et all ()2010D at low cluster masses, which pro- 
vides a reasonable description of the groups as well. This 
continues to agree well with the smaller groupings, all the 
way down to pairs of stars. 

While the small groupings follow the trend remarkably 
well, it is worth nothing that the range of values they can 
span becomes increasingly restricted at lower group sizes. 
The dashed line indicates a 1-1 relationship, the upper 
limit for any group size, representing all mass residing in 
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the most massive group member. For every number of 
group members, there is also a lower hmit, with a slope 
of the inverse of the number of members; this lower hmit 
hes parallel and below the 1-1 line on the plot. Even 
with this consideration of the smaller range possible for 
small groupings, the correlation between the maximum 
and total cluster mass is quite good, which suggests that 
properties of stellar groups may extend to systems with 
fewer than ten members. This could be interpreted as 
evidence that the masses of stars in the small groupings 
are consistent with a random sampling of the IMF or an 
IMF-like distribution, even though the total number of 
members is too small to measure this directly. Higher ac- 
curacy mass estimates would be needed to determine the 
consistency of the scatter in the data with that expected 
from pure random sampling. 

6. DISCUSSION 

In clustered regions, the tendency to find more massive 
stars in regions of higher surface density is linked to mass 
segregation, since higher surfa ce densities correspond to 
smaller cluster-centric radii. iHillenbrand fc HartmamJ 
(jT998) observed mass segregation in the ONC for stars 
with masses above 1-2 Mq at all radii, and some evidence 
for mass segregation at even smaller masses in the inner 
radii; a more recent analysis of the ONC data with a 
MST-based technique has found strong evidence of mass 
segregation only down to 16 M0, w ith weaker evidenc e 
of mass segreg ation down to 5 Mg ("All ison et all 120091 ). 
Simi larly, mass segregat ion has beeri observed for IC348 
(e.g. iMuench et alll2007|) . In iPaper^ . we analyzed mass 
segregation from the perspective of the groups, and found 
that the most massive group member, and in a few cases, 
also the second most massive group member, tended to 
be centrally located, reminiscent of mass segregation. 

In this work, we extend the analysis to include stars 
not in groups, and find that the differences in mass distri- 
butions apply not only to clusters, but also to dense and 
sparse groups. This tendency suggests a kind of "mass 
enrichment" associated with stellar surface density, but 
not "mass segregation" due to location within a group. 
We find that massive stars tend to be more prevalent in 
higher stellar surface density environments and within 
stellar groupings, well below the surface densities typi- 
cally associated with clusters. Perhaps the most surpris- 
ing result is that this trend is not only strong in IC348 
and the ONC where it had been previously observed, but 
is also seen in the relatively dispersed stellar population 
of Taurus. Given the tendency for more massive sources 
to be found in more clustered environments in both the 
clustered 1C348 and the dispersed Taurus regions, it is 
puzzling that a similar trend was not seen in Chal, whose 
environment appears to be intermediate between these 
two extremes. Two of the three groups identifi ed in Ch al 
did, however follow the general trend found it iPaper ll of 
having a centrally-located most massive group member. 

Stellar motion is the final piece of the puzzle not yet ex- 
amined. Since all of the regions are quite young, the stars 
should not have had a chance yet to move a significant 
distance away from their birth sites. We consider the 
available proper motion data for all of the early type stars 
in our sample to verify this assumption in Appendix [Xj 
Within the current observational uncertainties, there is 
no evidence for stellar motion to have influenced our re- 



sults. 

How do the differences in the stellar distributions in 
more versus less clustered environments arise? Are more 
low mass stars preferentially born in lower density envi- 
ronments, or are they former members of small groups 
that were cast out at an early age? Our data show sug- 
gestive signs that the excess in massive stars in more 
clustered environments may not be a uniform excess, but 
may be clustered around several spectral types, specifi- 
cally an excess of late B types and late K types, cor- 
responding to masses around 3 and IM0 at 1 Myr. If 
this observation were borne out in a larger dataset with 
higher statistical significance, that would suggest that a 
more complicated physical process related to formation 
mechanisms is responsible. It is difficult to see how a sim- 
ple picture of ejection, for example, would have such spe- 
cific preferred mass scales. If instead there is a roughly 
uniform excess of massive stars in more clustered en- 
vironment, this could be more easily attributed to the 
formation process. The fact that many stars form in a 
group suggests that the group environment has advanta- 
geous conditions relative to the field environment, such 
as a gas resevoir which allows for higher accretion rates 
or longer accretion timescales, which would be favourable 
for forming higher mass stars. 

The differences seen between the masses of clustered 
and isolated young stellar objects suggests that the initial 
mass function (IMF) is not independent of the physical 
conditions of star formation. Differences in the IMF may 
not always be so apparent, as the primary contributions 
are likely to be from the most 'successful' star-forming 
regions, i.e., more clustered environments. 

7. CONCLUSION 

Complementary to the analysis presented in IPaper J . 
we compare the more and less clustered populations of 
stars found in four nearby young star-forming regions - 
Taurus, Lupus3, Chal, and 1C348. The clustered young 
stellar population in the ONC is also used as a compari- 
son. We find the following results: 

1. Stars in the most clustered of regions, IC348 and 
the ONC, have statistically significant differences 
between the populations associated with high and 
low surface density and also with groups and iso- 
lation. The stars in high surface density envi- 
ronments or groups tend to be skewed overall to 
more massive members at the 98% or higher con- 
fidence level. These statistics are driven by by dif- 
ferences in the populations at the low-mass end, al- 
though differences are seen out to the earliest spec- 
tral types as well. The typical difference in masses 
between the high surface density or grouped stars 
and the low surface density or isolated stars is ~11- 
18% in IC348. 

2. Despite being the classic example of isolated, dis- 
tributed star formation, the stars in Taurus follow 
the same trends as the clustered regions. Stars in 
higher surface density environments in Taurus tend 
to have larger masses than the lower surface den- 
sity environment stars at the 99% confidence level, 
and stars in groups tend to have larger masses than 
those found in isolation at the 96% confidence level. 



11 



These statistics are again primarily driven by dif- 
ferences in the populations around spectral type M, 
although differences are seen at the earUest spectral 
types. The typical difference in masses between the 
high surface density or grouped stars and the low 
surface density or isolated stars is ~4-18%. 

3. The Lupus3 catalog is not large enough to allow 
for any strong conclusions to be made. 

4. Chal does not follow the trends seen in the other 
regions, and shows no evidence for more early type 
stars in g rouped o r higher surface density environ- 
ments. In iPaper ]| . two of the three groups in Chal 
did follow the general trend of a centrally-located 
most massive member, however. 

5. The combined mass distribution from all four re- 
gions is consistent to what is observed in the 
evolution-corrected field star mass function, or, 
similarly, the IMF. 

6. In all four regions, the relationship between the 
mass of the most massive group member and the 
total group mass seen in clusters and larger groups 
continues to arbitrarily small group sizes. 



7. Stellar motion does not appear to be responsible 
for the ejection of all the isolated early type stars, 
given the young ages of the star-forming regions, 
and the large separations of some of the early type 
stars from any stellar groupings. 

The complete stellar catalogues available for these four 
nearby star-forming regions allow the properties of stars 
in small groups and isolation to be quantified to an ex- 
tent not previously possible. Although theories of mas- 
sive and clustered star formation tend to focus on large, 
high density systems, our results suggest that small, rel- 
atively sparse, nearby regions offer an alternative regime 
in which to test the models. The deep, uniform com- 
pleteness levels and lack of contamination and source 
confusion allow for the characterization of grouped and 
isolated sources at a level not possible in more distant 
clusters. 
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APPENDIX 
PROPER MOTION 

The current YSO positions must be assumed to be similar to their birth locations in order to attribute the differences 
in stellar populations to the star formation process. Significant stellar motions are known to occur in some systems, 
however. In large clusters, the ejection of even relatively massive members is seen in numerical simulations (e.g., 
iVine fc Bo nncU 2003), and in observations (Goodman & Arcc 2004). Smaller, sparser systems better corresponding 
to our stellar groups are less well-studied. N-body simulations of very small groups (with N~5) including the effects 
of gas accretion and drag indicate that single stars up to one fifth the mass of the total group mass can be ejected 
early in the group's evolution ( Delgado-Donate, Clarke, & Bate 2003; see Kiseleva et al 1998 for the case of no gas). 
iBate. Bonnell. fc BromnJ ()2003[ ) furthermore argue that larger dynamical perturbations occur for both single stars and 
binaries in larger systems where larger-scale star-star interactions occur. Despite the young age and low stellar surface 
density in our star-forming regions, these simulations suggest that stellar motions should be examined to determine if 
they have could have effected the currently observed YSO positions. 

We examine the likelihood of significant stellar motion for the earlier type stars in our sample, both from a theoretical 
perspective of the necessary energetics (Appendix [Xjl) and using the available proper motion data (Appendix El 2). 
For ease of analysis, we divide the YSOs into our standard iV > 10 groups, and consider only the motion of the early 
type stars. The groups in general have good correspondence to high surface density environments, so the results should 
be generally applicable. Due to their fewer numbers, the migration of a single early type star would have a larger 
impact on the resultant distributions than a later spectral type star; observationally, proper motions are only available 
for the brightest stars. 

Theory 
Movement Out of Groups 

We first examine the likelihood that the isolated early type stars formed within their nearest group and subsequently 
migrated to their current locations. The escape speed required to leave each group is roughly 



2GMu 



R 



(Al) 



where Mtot is the current total mass of the group, and R is the separation from the group centre from which the star 
would start its migration. We consider for two cases for R - a, star forming at a 'typical' group separation, or Rrmsi 
the standard deviation of all current group member separations, and a star forming at the group outskirts, or Rmax, 
the current maximum group member separation. The values of R and Vesc for each group are given in Table 7. 

For each isolated early type star, we calculate the star's separation to both the centre and outskirts of its nearest 
group, using the separation to the closest group member as a proxy for the latter. For 1C348, separations only to 
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the main cluster are considered, since it strongly dominates the dynamics of the region; LupusS has no isolated early 
type stars, and is not considered. The separations are listed in Table 8; for stars 1 Myr old, each separation (in pc) 
is approximately equal to the velocity (in km s~^) required to move that distance. The ratios between the velocity 
required to reach the group centre or edge and the associated escape velocity are also given in Table 8, along with 
whether each star is more massive than all of the current group members. Figure [14] shows the separations of all 
isolated early type stars in Taurus from their closest group centre (top panel) and group outskirts (bottom panel) in 
right ascension and declination, with darker shades indicating earlier spectral types. The mean values of Rmax and 
Rrms for all of the groups in Taurus are shown in the bottom left corner of the top panel. The concentric annuli 
indicate separations of 0.5, 1, and 2 pc respectively. The other regions show similar behaviour, with some early type 
stars well-separated from the nearest group, and others located much closer. 

Quantitatively, for a plausible dynamical 'ejection' from the group, the isolated star should be less massive than the 
most massive member currently in the group, and the velocity required to take the star to its current location (assuming 
an age of 1 Myr) should not be substantially larger than the escape speed for the group. With these considerations, we 
see that all of the isolated early type stars in IC348 could have moved to their current locations, as well as few of the G 
stars in Taurus and Chal. Taurus-312, Taurus-352 and ChaI-209 and Cha-210, all appear too isolated and / or massive 
to be potential migration candidates. We stress that the migration candidate stars undoubtedly include a significant 
fraction which did not undergo migration. Our simple analysis uses only separations on the plane of the sky, while 
the 3D separation could be much larger. Even if the 3D separation is reasonable for migration, it does not mean that 
migration necessarily occured. Further constraints using available proper motion data are given in Appendix 1X12. 

Movement Into Groups 

For completeness, we also look for indications that the early type stars currently in groups migrated there after 
forming in isolation. Table 9 gives the separation between each sta r and th e group centre and edge, along with 
whether it is the most massive member of its group. As was found in iPaper^ . the early type stars tend to be found 
near the centre of their group, making it unlikely that they migrated there in only a ~1 Myr timescale. A few stars, 
such as Taurus-350 and Taurus-351 are the most likely candidates to have formed either outside the group, or on the 
outskirts of the group. 

Proper Motion Data 

We now use proper motion measurements in conjunction with the theoretical considerations above to constrain the 
likelihood of stellar migration. 

To fi nd the prop er motion of each s tar, w e searched the recent all-sky proper rn otion catal ogs of I Zacharias et all 
([20Toh . [R6ser et all (poTol ) . iRoser et al (pOOl . Fran Leeuwen I ([200l . IZacharias et all ([200l . and iKharehenko . L20m . 
In a non-negligible number of sources, the proper mo t ion li sted for the same source in multiple catalogs differed by 
more than the errors listed. In IC348. [Zacharias et all (|201Cf) list a subset of the stars as moving with a substantially 
different proper motion in RA (~25 mas yr""'^) than is given in the other catalogs (most give negative values). In 
order to minimize bias from discrepant catalogs, we use a two-step process to estimate the 'best' proper motions for 
each early type star, and caution that these values are still highly uncertain. We first calculate the median proper 
motion for each star, and then eliminate all measurements which differ by more than 15 mas/yr from this value (the 
quoted errors are typically much smaller than this). With the remaining measurements, we calculate the weighted 
mean (weighting by the error of each measure), which we will refer to as the 'best' value, and error. We also calculate 
the range spanned by all of the measures remaining after the first step; this range of 'good' proper motion measures 
provides a more realistic sense of the possible range in values than the formal error. Tables 10 and 11 list the 'best' 
proper motion, associated error, and range of possible values for the early type stars found in isolation and groups 
respectively. Note that the proper motion measurements are all given relative to each star's nearest or associated 
group. 

The proper motion of each early type star needs to be considered in relationship to the nearest or as sociated group. 
Ta ble 7 lists the proper motion we adopt for each group. In Taurus, we use the values given in Tubman et all ([2010); 
the iLuhman et all ()2010[ ) proper groups genera lly match well with the groups we identify, although both our Groups 
2 and 3 match the same ILuhman et all ( 20101) group. For Lupus3 and Chal, the group proper motions have been 
calculated in the same manner as for Taurus but are not yet published (E. Mamajeck, priv. comm.. Mar 9, 2011). No 
statistically significant differences were seen for the stars in the northern and southern parts of Chal, so the value for 
the entire population is used. A similar proper motion measurement has not yet been made for IC348, so we estimate 
this value using the weighted mean of the individual 'best' proper motions we calculate for group members; the error 
given is the standard error. 

Movement Out of Groups 

We estimate where each isolated star would have been 1 Myr ago relative to its nearest group, using the range of 
good proper motion values. The range of possible positions is often quite large; we list the smallest separation from the 
group in this range, SepiMyr, in Table 10 to indicate whether the star could possibly have originated from the group, 
given the available proper motion data. In IC348, the results are particularly uncertain, given both the larger scatter 
between measurements discussed earlier, as well as the small angular scale spanned by the stars (small relative motions 
have a larger impact on the final relative position). Given these caveats, we see that about half of the isolated early 
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type stars have proper motions consistent with having originated in a group, while the other half do not. Including the 
theoretical considerations in Appendix [^1, the stars which appear to be the most promising candidates for migration 
are Taurus-97, ChaI-91, and aU of the stars in IC348 barring IC348-7 and IC348-11. 

Movement Into Groups 

We similarly estimate where each grouped early type star would have been 1 Myr ago relative to its current group; 
these values are given in Table 11. Note that since SepiMyr measures the predicted minimum separation to the group 
centre 1 Myr ago, small, but non-zero values are consistent with the star being a group member in the past. With 
this consideration, the only stars whose proper motions suggest they formed outside the group and migrated in are 
Taurus-58, Taurus-351, Lupus3-22, Lupus-23, ChaI-154, and IC348-18. 

Upon closer consideration of the data, all of these sources appear to be consistent with originating in their current 
groups. Proper motion measurements consistent with grou p membersh i p are foun d for both Taurus-5 8 and Taurus- 
351 in catalogs which were excluded from the 'good' values (jRoser et~ail (|2008[ ) and iKharchenko I ()200in respectively); 
Taurus-351 is also a binary star, which may contribute to the discrepancy between catalogs. The range of good proper 
motion measures for both of the Lupus3 stars, ChaI-154, and IC348-18, are nearly compatible with the group's proper 
motion; given the scatter between catalog values, it seems likely that the errors are somewhat underestimated. 

Summary 

Significant stellar motion appears unlikely for most of the early type stars in our sample. A precise determination 
is difficult, given the large uncertainties in the proper motion measurements, the unknown line of sight distances, and 
precise ages of the stars. Isolated early type stars in IC348 have the highest likelihood of having originated in a group, 
since this region is the most compact and tends to have the most uncertain proper motion measurements. In Taurus 
and Chal, a few G-type stars currently found in isolation are the most likely candidates for having formed in a group. 
None of the stars currently in groups show compelling evidence for having formed in isolation and migrating inward. 
The sources which may belong to a different category than when they formed compose a small number of the total 
number of early type sources, and have primarily slightly later spectral types, where the total number of sources is 
larger. Our statistical analyses are therefore expected to remain similar even with proper motion considerations. 

Notably, the isolated FO and B6.5 stars in Chal appear unlikely to have originated in a group, suggesting that, while 
rare, some massive young stars may indeed form in isolation. 

EFFECT OF GROUP DEFINITION 

In Sectio n HI stars were classified as belonging to a group or being in isolation based on the stellar groups defined 
in iPapeF^ . i.e., stellar groups needed more than ten members, each separated from their nearest neighbour by less 
than the critical length Lcrit- We examine the effect of adopting a different minimum number of members to define 
'grouped' and 'isolated' YSOs in Appendix [Bll and[Bl2. and finally, using a different Lcrit in[B]3. 

Smaller Groups 

We re- run our comparisons of grouped and isolated YSO spectral types using minimum group sizes of N >10, 8, 6, 
4, and 2. We tend to find similar statistical results for the smaller minimum group sizes, with the variation increasing 
at the smallest group sizes. Table 6 gives the KS2 and MW probabilties for a minimum group size of iV > 2. The 
grouped and isolated stars in ONC are clearly still distinct with this drastic change in group defintion. In Lupus3 and 
Chal, as before, the probabilities do not allow us to make definitive conclusions. In IC348, the grouped and isolated 
populations appear to be much more similar with N > 2 than N > 10, which appears to be due to the small number 
of stars which are classified as isolated (the grouped stars are nearly seven time more numerous than the isolated stars 
in this instance). For > 6, the KS2 and MW statistics are in better agreement with the standard A^ > 10 case. 
In Taurus, the populations of grouped and isolated stars actually appear to be more distinct for N > 2 than with 
N > 10. This is in agreement with the finding in Section 3 that small groupings of stars in Taurus are distinct from 
the truly isolated stars. 

Larger Groups 

Conversely, we examine the effect of increasing the minimum group size. Since there are few groups of any given 
size, we examine minimum group sizes of A^ > 15, 20, 25, 40, and 55, chosen to match the discrete group size values in 
our sample. These values also encompass group sizes in which the number of grouped and isolated sources are roughly 
equal, minimizing the statistical uncertainty in comparisons. The grouped and isolated source counts are closest for 
10 < A^ < 16 (Taurus), 12 < A < 36 (Chal), 11 < A^ < 186 (IC348), 37 < A^ < 410 (ONC-c), while in ONC-f, A^ < 10 
would be required. Lupus3 has only one group with A > 10, so cannot be tested for an increase in the minimum group 
size. 

Table 6 shows the KS2 and MW probabilities for A^ > 15. As the minimum group size, A^ is increased, the KS2 
probabilities that the 'grouped' and 'isolated' spectral types are drawn from the same parent distribution tend to 
increase slowly with small increases in and then rapidly increase once a significant fraction of the original A^ > 10 
grouped stars become classified as 'isolated'. This transition tends to occur at or slightly below the value of A^ which 
makes the number of sources classified as grouped or isolated roughly equal. From this, we conclude that a slight 
increase of the minimum group size does not have a significant impact on our results. Increasing the minimum group 
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size significantly (~20 or higher, depending on the region), however, serves to wash out the differences seen between 
grouped and isolated sources. 

Different Lcrit 

Lcrit, the maximum separation between group m embers, is the second variable which controls whether sources are 
classified as belonging t o a group or not. In lPaper"^. we fou nd an uncertainty of < 10% in Lcrit, using the definition in 
iGutermuth et all (|2009( ) for Lcrit- Appendix D.l of jPaper II examined the effect of a varying Lcrit on group membership, 
and found typically only a handful of late type stars would be included or excluded from all of the groups in a region 
with a different Lcrit- Changes much larger than the measurement error (at least 40-60%, but usually at least double 
the measured value) are required to cause significant changes to the classification of sources. We therefore conclude 
that the choice of Lcrit does not have a strong influence on the resulting spectral type distributions or the statistical 
measures that we present. 
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Table 1 

Statistics for YSOs separated by surfa<;e density 



n = 4 n = 9 



Region 


(%) 


(%) 


Nearly 

max a cutoff 


Median S'* 
low up 


KS2=' 
(%) 


(%) 


max (T 


cutoff 


Median S^* 
low up 


Taurus 


5 


99 


1.4 


GO,AO 


0.65 


44.0 


54 


92 


1.4 


GO,AO 


0.61 


15.8 


LupusS 


96 


75 


1.4 


G9-G0 


1.80 


154 


96 


62 


1.4 


G9-G0 


2.56 


125 


Chal 


25 


32 


0.7 


A5,A0 


4.95 


140 


25 


25 


1.5 


G9 


4.97 


1.38 


IC348 


3 


99 


3.5 


GO 


48.3 


667 


0.3 


98 


1.8 


GO 


48.2 


491 


ONCl-c° 


0.1 


99.97 


4.9 


B5 


28.0 


379 


0.005 


99.995 


3.0 


B5 


27.7 


276 


ONCl-f= 


0.8 


98 


7.1 


B5 


38.7 


977 


5 


90 


7.8 


B5 


37.4 


910 


Combined^ 


0.6 


99.8 


3.5 


GO 






5 


97 


2.4 


GO 






T+L+IS 


0.2 


99.98 


4.5 


GO 






1 


99 


2.4 


GO 






Combined global*^ 


68 


55 


4.1 


GO 


2.72 


285 


0.08 


0.1 


3.5 


GO 


2.24 


234 


T+L+I global' 


28 


57 


5.3 


GO 


2.24 


345 


6 


29 


4.1 


GO 


2.02 


292 


L+C+I global^ 


9 


98 


6.5 


GO 


10.9 


425 


22 


91 


6.0 


GO 


11.3 


339 



^ The two-sample Kolmogorov-Smirnov probability that the spectral types of stars in high and low surfa<;e density environments 
are drawn from the same parent distribution. 

^ The Mann- Whitney probabiUty that the spectral types of stars in high surface density environments tend to be earUer than 
the spectral types of stars in low surface density environments. 

The maximum (sigma) excess of early type stars in high versus low surface density environments, and the early spectral type 
cutoff at which this occurs. 

The median surface density, in pc~^, of YSOs in each of the two surface density categories. 
" The ONC catalog using a conservative cut (-c) and the full catalog with no cuts (-f). 
^ Statatistics for the combined Taurus, Lupus3, Chal, and IC348 catalogs. Sec Section 3.4 for details, 
s Statistics for the combined Taurus, Lupus3, and 1C348 catalogs. See Section 3.4 for details. 

^ Statistics for Taurus, Lupus3, Chal, and IC348, separated by the global surface densities. See Section 3.4 for details. 
' Statistics for Taurus, Lupus3, and IC348, separated by the global surface densities. See Section 3.4 for details. 
J Statistics for Lupus3, Chal, and IC348, separated by the global surface densities. See Section 3.4 for details. 



Table 2 

Effect of Spectral Type Uncertainty on Surface Density Statistics 



KS2 Statistic 







n = 4 






n = 9 




Region 


median " 


90th percentile " 


median " 


90th percentile " 




(%) 


(%) 


(%) 


{%) 


(%) 


(%) 


Taurus 


8 


4 


11 


54 


34 


76 


Lupus3 


96 


46 


96 


96 


74 


96 


Chal 


60 


25 


96 


60 


25 


96 


IC348 


3 


0.8 


8 


0.5 


0.1 


2 


ONCl-c'' 


0.05 


0.02 


0.1 


0.003 


0.001 


0.007 


ONCl-f'' 


0.6 


0.3 


1 


3 


2 


6 


Combined'^ 


0.7 


0.2 


2 


6 


2 


14 


T-hL+ld 


0.2 


0.05 


0.7 


1 


0.5 


4 


MW Statistic 


Taurus 


99 


99 


99.6 


92 


89 


94 


Lupus3 


75 


60 


96 


62 


48 


74 


Chal 


35 


26 


45 


27 


19 


36 


1C348 


99 


98 


99.5 


98 


96 


99 


ONCl-c'" 


99.96 


99.9 


99.98 


99.99 


99.98 


99.997 


ONCl-f*^ 


97 


95 


98 


88 


83 


92 


Combined'^ 


99.8 


99.6 


99.9 


98 


96 


99 


T-I-L-I-I'l 


99.97 


99.9 


99.99 


99 


99 


99.7 



a The median and the range that 90% of the values he within for the KS2 or MW statistic found with the 
random addition of an error of half a spectral subtype. 

^ The ONC catalog using a conservative cut (-c) and the full catalog with no cuts (-f). 
" Statatistics for the combined Taurus, LupusS, Chal, and IC348 catalogs. 
Statistics for the combined Taurus, LupusS, and ICS48 catalogs. 
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Table 3 

Statistics for Grouped and Isolated YSOs 



Region 


KS2'' 












(%) 


(%) 


max a 


cutoff 


Taurus 


32 


96 


1.4 




GO 


LupusS 


95 


38 


1.4 




G9-G0 


Chal 


11 


12 


-0.4 




GO 


IC348 


0.8 


98 


1.2 




G9 


ONCl-cd 


0.02 


99 


1.2 




B5 


ONCl-Pi 


6 


86 


2.0 




B5 


Combined'* 


5 


94 


1.5 




G9 


T+L+lf 


1 


99 


1.8 




GO 



^ The two-sample Kolmogorov-Smirnov probability that the spectral types of grouped and isolated stars are 
drawn from the same parent distribution. 

^ The Mann- Whitney probability that the speetral types stars in groups tend to be earlier than the spectral 
typos of stars in isolation. 

The maximum (sigma) excess of early type stars in groiips veruss isolation, and the early spectral type 
cutoff at which this occurs. 

A comparison to the ONC catalog, using a conservative cut (-c) and the full catalog with no cuts (-f). 
" Statistics for the combined Taurus, LupusS, Chal, and IC348 catalogs. 
' Statistics for the combined Taurus, LupusS, and IC348 catalogs. 



Table 4 

Effect of Spectral Type Uncertainty on MST Group Statistics 







KS2 Statistic 






MW Statistic 




Region 


median " 


90tfi percentile " 


median " 


90th percentile 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


Taurus 


28 


14 


47 


96 


94 


97 


Lupus3 


89 


48 


99.7 


38 


26 


53 


Chal 


28 


9 


63 


15 


9 


22 


IC348 


0.2 


0.04 


0.8 


98 


96 


99 


ONCl-c'' 


0.2 


0.04 


0.9 


99 


99 


99.6 


ONCl-f'' 


5 


3 


9 


85 


79 


90 


Combined"^ 


8 


4 


19 


95 


92 


97 


T-l-L-l-I'' 


2 


0.6 


5 


99 


99 


99.7 



The median and the range that 90% of the values lie within for the KS2 or MW statistic found with the 
random addition of an error of half a spectral subtype. 

^ The ONC catalog using a conservative cut (-c) and the full catalog with no cuts (-f). 
Statatistics for the combined Taurus, LupusS, Chal, and IC348 catalogs. 
Statistics for the combined Taurus, LupusS, and ICS48 catalogs. 



Table 5 

Statistics for Individual Groups in Taurus and Chal 



Region 


Group 


KS2^ 








(%) 


(%) 


Taurus 


1 


60 


64 


Taurus 


2 


88 


67 


Taurus 


3 


71 


48 


Taurus 


4 


3 


99.8 


Taurus 


5 


8 


89 


Taurus 


6 


99 


67 


Taurus 


7 


73 


83 


Taurus 


8 


54 


61 


Chal 


1 


42 


10 


Chal 


2 


12 


5 


Chal 


3 


46 


72 



The two-sample Kolmogorov-Smirnov probability that the spectral types of grouped and isolated stars are 
drawn from the same parent distribution. 

^ The Mann- Whitney probability that the spectral types stars in groups tend to be earlier than the spectral 
types of stars in isolation. 
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Table 6 

Statistics for Groupings and Isolated YSOs with N > 2 and > 15 



Region 


KS2=- 
(%) 


TV > 2 


(%) 


KS2'' 
(%) 


TV > 15 


MW'' 
(%) 


Taurus 


7 




95 


42 




90 


LupusS 


37 




13 


95 




38 


Chal 


19 




68 


23 




26 


IC348 


17 




63 


0.3 




98 


ONCl-c= 


0.8 




91 


0.8 




98 


ONCl-F^ 


0.02 




81 


0.4 




98 



^ The two-sample Kolmogorov-Smirnov probability that the spectral types of stars in any size of grouping 
and very isolated stars are drawn from the same parent distribution. 

The Mann- Whitney probability that the spectral types of stars in any size of grouping tend to be earlier 
than the spectral types of very isolated stars. 

" A comparison to the ONC catalog, using a conservative cut (-c) and the full catalog with no cuts (-f). 



Table 7 

Group Properties 



Region 


Grp 


Rrms^ 


Rmax'^ 


Vesc(Rrms)'' 


Vesc(Rmax)'' 


Ma'' 






Refs <^ 




# 


(pc) 


(pc) 


(km s^^) 


(km s-i) 


(mas/yr) 


(mas/yr) 


(mas/yr) 




Taurus 


1 


0.24 


0.82 


0.62 


0.33 


6.9 


-22.3 


1 


LIO 


Taurus 


2 


0.29 


1.13 


0.68 


0.34 


6.0 


-26.8 


1 


LIO 


Taurus 


3 


0.40 


1.31 


0.42 


0.23 


6.0 


-26.8 


1 


LIO 


Taurus 


4 


0.35 


1.38 


0.75 


0.38 


10.0 


-17.6 


1 


LIO 


Taurus 


5 


0.23 


0.69 


0.56 


0.32 


5.5 


-21.9 


1 


LIO 


Taurus 


6 


0.42 


1.64 


0.60 


0.30 


6.7 


-17.7 


2 


LIO 


Taurus 


7 


0.48 


1.85 


0.54 


0.27 


4.5 


-21.3 


2 


LIO 


Taurus 


8 


0.26 


0.88 


0.64 


0.35 


3.9 


-23.4 


1 


LIO 


Lupus3 


1 


0.24 


1.06 


0.80 


0.38 


-9.4 


-21 


1.1 


EMI 


Chal 


1 


0.16 


0.48 


0.44 


0.26 


-20.9 


1.3 


0.7 


EM2 


Chal 


2 


0.21 


1.01 


1.27 


0.59 


-20.9 


1.3 


0.7 


EM2 


Chal 


3 


0.17 


0.61 


1.04 


0.55 


-20.9 


1.3 


0.7 


EM2 


IC348 


1 


0.13 


0.58 


2.73 


1.29 


5.6 


-6.6 


0.8 


KM 


IC348 


2 


0.04 


0.15 


0.79 


0.43 


5.6 


-6.6 


0.8 


KM 



The rms and maximum separation of group members from the group centre. 
^ Escape velocity for the group, calculated at both Rrms and Rmax 
Proper motion measured for each group (RA, dec, and approximate standard error) 

References for proper motion measurement. LIO: ILuhman et all 12010) : Groups 2 and 3 both lie within Luhman's proper 
motion Group II. EMI: E. Mamajek, priv. comm, Mar 9, 2011. EM2: E. Mamajek, priv. comm. Mar 9, 2011; no statistically 
significant difference is seen between the northern and southern stars, so the average value for the entire region is adopted. KM: 
Calculated by the authors; more details in text. 
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Table 8 

Required Motions for Isolated YSOs 



IXC^lOll 


Src 


Spec 




OtJpcLl txtlUli 


oc^ai atiuii 


V / Vesc'' 


V / Vesc'^ 


M < 




M- 

W 


lype 




CcntrG (pc) 




Centre 


riiCige 


M ? 


Taurus 


97 


G5 


2 


1.69 


0.84 


2.53 


2.49 


Y 


Taurus 


160 


G8 


6 


2.50 


0.95 


4.26 


3.20 


Y 


Taurus 


233 


G8 


4 


1.35 


0.86 


1.83 


2.34 


Y 


Taurus 


312 


G5 


4 


9.77 


8.51 


13.29 


23.16 


Y 


Taurus 


352 


A2 


8 


2.19 


1.67 


3.49 


4.92 


Y 


Chal 


91 


G9 


3 


0.63 


0.23 


0.62 


0.43 


Y 


Chal 


209 


FO 


1 


2.56 


2.21 


5.95 


8.87 


N 


Chal 


210 


B6.5 


3 


1.38 


0.92 


1.35 


1.70 


N 


IC348 


3 


AO 




0.94 


0.44 


0.35 


0.35 


Y 


IC348 


7 


AO 




0.49 


0.18 


0.18 


0.14 


Y 


IC348 


8 


A2 




0.48 


0.17 


0.18 


0.14 


Y 


IC348 


11 


G4 




0.76 


0.50 


0.29 


0.39 


Y 


IC348 


17 


A4 




0.93 


0.44 


0.35 


0.35 


Y 


IC348 


19 


Gl 




0.64 


0.37 


0.24 


0.29 


Y 


IC348 


21 


G5 




0.87 


0.50 


0.32 


0.39 


Y 


IC348 


25 


A4 




0.62 


0.37 


0.23 


0.29 


Y 



For IC348, all comparisons arc to IC348-niain; sources 3, 7, and 21 lie closer to IC348-Nortli. 
^ Ratio of the velocity required for the source to move to its current location from the group centre to the 
escape velocity at Rrms. 

Ratio of the velocity required for the source to move to its current location from the nearest current group 
member to the escape velocity at the group edge. 

Is the source less massive than the most massive group member? 



Table 9 

Required Motions for Grouped YSOs 



Region 


Src 


Spec. 


Group^ 


Separation 


Separation 


M = 




# 


Type 




Edge (pc) 


Centre (pc) 




Taurus 


58 


B9 


2 


1.01 


0.12 


Y 


Taurus 


255 


GO 


6 


1.27 


0.36 


Y 


Taurus 


287 


G5 


7 


1.44 


0.40 


Y 


Taurus 


328 


B9 


8 


0.68 


0.20 


Y 


Taurus 


336 


G2 


8 


0.62 


0.26 


N 


Taurus 


350 


B9.5 


4 


0.43 


0.95 


N 


Taurus 


351 


B9 


4 


0.43 


0.95 


Y 


Lupus3 


21 


A7 


1 


1.03 


0.03 


N 


Lupus3 


22 


A3 


1 


1.02 


0.03 


Y 


Chal 


11 


G5 


2 


0.43 


0.58 


N 


Chal 


53 


B9 


3 


0.47 


0.13 


Y 


Chal 


116 


G2 


2 


0.63 


0.37 


N 


Chal 


121 


B9.5 


2 


0.63 


0.38 


N 


Chal 


154 


G8 


2 


0.27 


0.74 


N 


Chal 


182 


G7 


2 


0.87 


0.13 


N 


IC348 


1 


B5 




0.51 


0.07 


Y 


IC348 


2 


A2 




0.48 


0.10 


N 


IC348 


4 


FO 




0.34 


0.23 


N 


IC348 


5 


G8 




0.16 


0.41 


N 


IC348 


6 


G3 




0.37 


0.20 


N 


IC348 


9 


G8 




0.47 


0.10 


N 


IC348 


10 


F2 




0.37 


0.20 


N 


IC348 


12 


GO 




0.33 


0.24 


N 


IC348 


13 


A3 




0.33 


0.24 


N 


IC348 


16 


G6 




0.54 


0.04 


N 


IC348 


18 


A2 




0.48 


0.10 


N 


IC348 


28 


FO 




0.30 


0.27 


N 


IC348 


29 


Gl 




0.12 


0.45 


N 


IC348 


35 


GO 




0.32 


0.25 


N 


IC348 


359 


B5 




0.51 


0.07 


Y 


IC348 


360 


G8 




0.47 


0.10 


N 



For IC348, all comparisons are to IC348-main; sources 3, 7, and 21 lie closer to IC348-North. 
^ Is the source the most massive group member? 



Table 10 

Observed Relative Motions for Isolated YSOs 



Region 


Src 


,j a 
a. 




M(5 




Min fia^ 


Max /J-ct^^ 


Mill fis^ 


Max fi^^ 


1) 1 A/Tm. 

'-'^ tr 1 Ivl y 


Taurus 


97 


1.63 


1.14 


-1.05 


1.13 


-1.08 


4.02 


-3.90 


1.15 





Taurus 


160 


-5.70 


4.60 


-27.33 


4.58 


-32.19 


2.12 


-44.29 


-21.18 


14.6 


Taurus 


233 


-13.26 


2.42 


9.61 


2.42 


-15.68 


-10.84 


7.19 


12.03 


9.6 


Taurus 


312 


2.66 


1.19 


-1.70 


1.19 


-0.11 


4.70 


-3.29 


0.04 


5.9 


Taurus 


352 


1.25 


1.09 


-1.53 


1.10 


-0.97 


3.09 


-4.29 


0.28 





Chal 


91 


-0.40 


1.11 


1.11 


1.04 


-3.87 


3.11 


-2.81 


5.97 





Chal 


209 


-5.78 


0.80 


-11.47 


0.72 


-7.19 


-4.53 


-13.29 


-8.79 


6.4 


Chal 


210 


-1.90 


0.73 


0.12 


0.67 


-3.50 


-0.63 


-0.91 


1.30 


0.2 


IC348 


3 


-2.30 


0.93 


0.48 


0.68 


-3.94 


3.94 


-2.09 


1.41 





IC348 


7 


2.78 


1.16 


-1.93 


1.00 


0.76 


5.55 


-5.31 


-0.01 


1.6 


IC348 


8 


-0.04 


1.16 


-0.69 


0.98 


-1.63 


2.38 


-5.42 


1.23 





IC348 


11 


-8.38 


4.84 


-4.04 


4.80 


-22.45 


-2.74 


-18.90 


1.72 


4.6 


IC348 


17 


-2.54 


5.45 


-2.10 


4.41 


-30.92 


3.76 


-10.62 


34.51 





IC348 


19 


-5.34 


2.74 


-0.66 


4.70 


-9.10 


2.26 


-8.32 


15.31 





IC348 


21 


-5.49 


2.27 


2.90 


4.04 


-8.14 


1.16 


-7.02 


20.05 





IC348 


25 


-1.45 


3.83 


-0.16 


4.01 


-6.76 


4.06 


-8.82 


9.62 






The estimated 'best' proper motion of each source relative to its nearest group and the associated errors, both in mas yr~^. 
The range of 'good' proper motion measurements of each source relative to its nearest group, all in mas yr~^. 
The minimum separation to the nearest group centre 1 Myr ago (in pc), based on the range of relative proper motion 
measurements. 



Table 11 

Observed Relative Motions for Grouped YSOs 



Region 


Src 








Aw r'' 


Min ^la'° 


Max /ia*" 


Min ^ls'° 


Max fMs^ 


SepiMj/i 




# 


(mas yr ^) 


(mas yr ^) 


(mas yr ^) 


(mas yr ^) 


(pc) 


Taurus 


58 


-27.15 


8.12 


-2.02 


4.26 


-39.47 


-15.06 


-8.10 


9.73 


10.1 


Taurus 


255 


4.71 


2.21 


3.73 


2.21 


1.49 


7.56 


-0.22 


6.90 


0.6 


Taurus 


287 


-2.90 


6.23 


1.60 


6.23 


-9.13 


3.33 


-4.63 


7.83 





Taurus 


328 


-1.68 


1.13 


-1.17 


1.06 


-3.35 


1.24 


-2.93 


0.25 


0.03 


Taurus 


336 


0.79 


1.26 


-0.25 


1.19 


-5.78 


3.22 


-3.75 


3.04 





Taurus 


350 


-2.95 


1.44 


-0.64 


1.11 


-6.76 


2.33 


-3.36 


1.87 





Taurus 


351 


15.20 


2.09 


-1.66 


1.49 


12.41 


17.99 


-6.41 


2.05 


7.7 


LupusS 


21 


0.18 


1.17 


-2.91 


1.17 


-1.69 


4.44 


-7.43 


-0.96 


0.9 


Lupus3 


22 


3.10 


1.18 


-1.22 


1.16 


0.91 


5.42 


-3.25 


0.25 


0.9 


Chal 


11 


2.26 


1.01 


-1.79 


1.78 


0.51 


7.48 


-4.75 


1.18 


0.8 


Chal 


53 


-0.54 


0.84 


-1.91 


0.76 


-1.99 


1.46 


-3.29 


0.48 





Chal 


116 


-3.29 


1.36 


3.06 


1.28 


-5.65 


0.38 


-1.04 


8.05 





Chal 


121 


1.42 


0.78 


0.41 


0.71 


-0.75 


3.25 


-1.37 


1.77 





Chal 


154 


-0.69 


1.20 


3.52 


1.17 


-1.89 


0.51 


2.35 


4.69 


1.4 


Chal 


182 


-0.44 


2.73 


6.66 


2.38 


-10.15 


8.87 


-11.21 


16.71 





IC348 


1 


-1.03 


1.15 


-0.15 


0.97 


-8.04 


0.64 


-6.77 


16.71 





IC348 


2 


1.39 


1.44 


0.23 


1.32 


-2.02 


3.95 


-4.46 


2.56 





IC348 


4 


2.52 


1.21 


-0.24 


0.87 


-4.74 


5.30 


-2.31 


2.56 





IC348 


5 


-1.25 


3.73 


-2.51 


3.60 


-8.96 


6.26 


-13.88 


6.92 





IC348 


6 


-6.58 


3.55 


-2.21 


4.59 


-12.77 


1.86 


-9.72 


22.11 





IC348 


9 


6.29 


1.59 


3.78 


3.08 


-10.15 


8.31 


-11.52 


9.34 





IC348 


10 


0.92 


1.10 


3.41 


1.12 


-1.84 


3.80 


-0.43 


5.09 





IC348 


12 


-4.80 


5.66 


-2.13 


4.94 


-10.46 


0.86 


-14.11 


4.02 





IC348 


13 


-4.80 


5.66 


-2.13 


4.94 


-10.46 


0.86 


-14.11 


4.02 





IC348 


16 


-3.20 


3.10 


-3.30 


2.94 


-6.30 


-0.10 


-6.24 


-0.36 


0.5 


IC348 


18 


2.91 


1.50 


-1.12 


1.36 


0.93 


4.90 


-3.00 


0.77 


1.4 


IC348 


28 


0.37 


1.46 


0.85 


1.40 


-5.85 


4.10 


-2.35 


2.57 





IC348 


29 


-5.20 


6.95 


-6.00 


6.92 


-12.15 


1.75 


-12.92 


0.92 





IC348 


35 


-3.83 


3.45 


-1.65 


3.65 


-10.26 


0.88 


-6.63 


26.31 





IC348 


359 


-1.01 


1.10 


-0.25 


0.92 


-12.81 


0.64 


-3.23 


16.71 





IC.318 


360 


6.29 


1.59 


3.78 


3.08 


-10.15 


8.31 


-11.52 


9.34 






* The estimated 'best' proper motion of each source relative to its group and the associated errors. 
The range of 'good' proper motion measurements of each source relative to its group. 
The minimum separation to the nearest group centre 1 Myr ago. 
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Figure 1. The fraction of stars with spectral types earlier than each value for stars in the highest (blue) and lowest (orange) third of 
stellar surface densities, computed using n = 4. The vertical lines indicate the Poisson error at each spectral type, and the numbers in 
braces in the legend indicate the total number of stars in each category. The inset plot shows the cumulative number in log space. Stars 
in Taurus and LupusS are shown in this figure. 



21 

Choi 

I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I— 




LO MO KO GO FO AO BO 

Spectrol type 

IC348 

—I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




LO MO KO GO FO AO 

Spectrol type 



Figure 2. The distribution of spectral types in the highest and lowest third of source surface densities in Chal and IC348. See Figure[T] 
for more details. 
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Figure 3. The distribution of spectral types in the highest and lowest third of source surface densities in the ONC. See Figure[T]for more 
details. 
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Figure 4. The distribution of spectral types in the highest and lowest third of surface densities combined over Taurus, LupusS, Chal, 
and IC348. The top panel shows the distributions when combined using the divisions in Figures [T] and [2] The bottom panel shows the 
distributions when divided using the highest and lowest third of the surface densities compared across all four regions. 
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Figure 5. The distribution of spectral types for grouped and isolated stars. See Figure [T] for the plotting conventions used. Stars in 
Taurus (top) and LupusS (bottom) are shown in this figure. 
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Figure 6. The distribution of spectral types for grouped and isolated stars. See Figure[T]for the plotting conventions. Stars in Chal (top) 
and IC348 (bottom) are shown in this figure. 
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Figure 7. The distribution of spectral types for grouped and isolated stars. See Figure [T] for the plotting conventions. Stars in the ONC 
are shown in this figure. 
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Figure 9. The distribution of spectral types in the two Taurus groups showing the most significant differences with the isolated populations 
(based on KS2 and MW statistics). The plot conventions follow those in Figure [T] 
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Figure 10. The cumulative distribution of mass for all four nearby star-forming regions combined. The top panel shows the high and 
low stellar surface density stars, while the bottom panel shows the grouped and isolated stars. Error bars indicate the Poisson errors. Also 
shown are predictions from MIO; for stars in grouped and isolated environments, scaled to a total fraction of 15% at the lower mass limit 
of O.SMq (black dashed and dotted lines respectively). The range of values spann ed by the analytic IKifl model are shown by the green 
asterisk (high density model) and red line (range of low density models). The lKlOl models are not directly comparable to the observations, 
however, the range spanned by the models should encompass the observations. 
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Figure 11. The cumulative mass distribution for the stars in ONC-c, for the high ve rsus l ow stella r surface density source s (t op panel), 
and grouped versus isolated sources (bottom panel) . Comparisons with predictions by IMIOI and IKICI are shown as in Figure 1101 note that 
the plot ranges differ from the previous figure. 
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Figure 12. The diflferential mass distribution for all four nearby star-forming regions. The top panel shows the high and low stellar 
surface density stars, while the bott om panel sh ows the grouped and isolated stars. Error bars indicate Poisson errors. Also shown is the 
field star mass function measured bv lReid et all | [2002f ) . corrected for stellar evolution (red line). The errors estimated bv lReid et all I I2002I) 
are shown offset from the bin centre for clarity. 
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Figure 13. The mass of the maximu m mass member in each stellar grouping versus the total group mass. The data and approximate 
low-mass slope in lWeidner et a] l|2010l ) are shown by the grey diamonds and dotted line respectively. The coloured circles show our results 
for groupings with six or more members (top) and two or more members (bottom). The dashed line shows a 1-1 relationship (the upper 
limit) . 
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Figure 14. The separation between each isolated early type star and its nearest group's centre (top) and the closest group member 
(bottom) in Taurus. Colours denote different spectral types, and the circles indicate separations of varying amounts. The thick horizontal 
lines show the mean Rrms and Rmax values for the groups. 



